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We search for direct CP violation in charmless hadronic B decays observed in a sample of about 
22.7 million BB pairs collected with the BABAR detector at the PEP-II asymmetric-energy e^e~ 
collider. We measure the following charge asymmetries: Acp^B^ r/' K^) = —0.11 ± 0.11 ± 0.02, 
Acp{B^ u}TT^) = -0.01 1 3'^ ± 0.03, Acp{B^ <j)K^) = -0.05 ± 0.20 ± 0.03, Acp{B^ ^ 

= -0.43 t aao ± 0-06, and Acp{ b'" ^ </> ''K'") = 0.00 ± 0.27 ± 0.03. 
PACS numbers: 13.25.Hw, ll.30.Er, 14.40.Nd 



The phenomenon of Charge-Parity (CP) symmetry vi- 
olation has played an important role in understanding 
fundamental physics since its initial discovery in the K 
meson system in 1964 [Q]. Soon after, it was recognized 



that the violation of CP symmetry was one of the fun- 
damental requirements to produce a matter-dominated 
Universe A significant CP- violating asymmetry in 
decays of neutral B mesons to final states containing 
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charmonium, due to interference between — mix- 
ing and direct decay amplitudes, has recently been ob- 
served 1^. As it has now been established Q that the 
CP-violating decays of the meson to tttt final states 
are due to CP violation in decay amplitudes as well as 
to K'^ — mixing, it is topical to search for "direct" 
CP asymmetries in B decays, which involve only direct 
decay amplitudes. These asymmetries are anticipated to 
be much larger in B decays than in K decays |^ . Direct 
CP violation would be measured as an asymmetry of B 
decay rates: 



A, 



CP 



TjB^ f)~TiB^ f) 

r{B ^ /) + r(p ^ /) 



(1) 



Charmless B meson decays are particularly interesting 
processes to search for direct CP violation because of 
the possible involvement of penguin (P) and tree (T) 
amplitudes of comparable magnitude. Substantial CP 
violation can thus arise in the Standard Model through 
interference of these terms H: 



Acp = 



2 |P| \T\ sinAcj) sinA(5 



P|2 + |r|2 + 2 |P| \T\ cos cos A(5 ' 



(2) 



where A(f) and AS are the differences in weak and strong 
phases. Because of the weak phase difference between 
the tree and penguin amplitudes, Acp is sensitive to 
the Cabibbo-Kobayashi-Maskawa matrix phases 7 = 
argf-KdK^/KA;] and a ^ ^rg[-V,,V;j V^,V:,]. 
The difference between the h ^ u tree and 6 — > s (& — > d) 
penguin amplitude weak phases A0 is 7 (a), as in the 
case of the decays B ttK, rj'K [B tttt, uitt). How- 
ever, large uncertainties in the strong phases, which can 
be calculated by certain models, weakens the quantita- 
tive relationship to the weak phases. Recent calculations 
based on effective theory and factorization predict asym- 
metries as large as ^10% Q. 

The measurement of direct CP violation in pure pen- 
guin modes, such as P ^ (pK^*\ is more sensitive to 
non-Standard-Model physics. In the Standard Model, 
the lack of a tree-level contribution results in an ex- 
pected Acp of no more than ~1% However, new 
particles in loops, such as charged Higgs boson or SUSY 
particles, would provide additional amplitudes with dif- 
ferent phases. Depending on the model parameters, Acp 
can be 30% or larger in such scenarios Complemen- 
tary searches for new physics would involve measure- 
ments of the time-dependent asymmetries in P ° de- 
cays to CP eigenstates, such as rj'Kg^j^y (j)Kg^j^y and 

(/) K^tt'^). Comparison of the value of sin 2/3 ob- 

tained from these modes with that from charmonium 
modes |^ can probe for new physics participating in pen- 
guin loops. In these measurements, direct CP violation 
in the decay becomes highly relevant and can be studied 
in the self-tagging modes discussed below. 



The CLEO experiment has reported a search for di- 
rect CP violation in B meson decays to ttK, rj'K , and 
ujTT 1^. In this paper we improve the precision of the 
measurements in the f]' K and wtt modes and extend 
the search to new modes. Measurements from BABAR 
of the B ttK charge asymmetries are presented else- 
where |ll|. Here we present measurements of the 
charge asymmetries in the following charmless B decays, 
for which branching fractions have been previously re- 
ported |l|: B^ ^ v'K^, B^ cj7r±, P± ^if^, 

P± (t>K*^, and *P'° (t>K^*°. The B flavor is de- 
termined by its charge, except for the (f> K *^ final state 
where the flavor is determined from the charge of the 

kaon from the K ^ K^n^ decay. 

The data were collected with the BABAR detector Q 
at the PEP-II asymmetric-energy e+e~ collider lo- 
cated at the Stanford Linear Accelerator Center. The 
results presented in this paper are based on data taken 
in the 1999-2000 run comprising an integrated luminosity 
of 20.7 fb^^, corresponding to 22.7 million BB pairs, at 
the T{4:S) resonance ( "on-resonance" ) and 2.6 fb~^ ap- 
proximately 40 MeV below this energy ( "off- resonance" ) . 
The T{4S) resonance occurs at the e+e~ center-of-mass 
(cm.) energy, -y/s, of 10.58 GeV. 

Charged particles are tracked and their momenta mea- 
sured with a combination of a silicon vertex tracker 
(SVT) consisting of five double-sided detectors and a 40- 
layer central drift chamber (DCH), both operating in a 
1.5 T solenoidal magnetic field. With the SVT, a posi- 
tion resolution near the interaction point of about 40 /im 
is achieved for the highest momentum charged particles, 
allowing the precise determination of decay vertices. The 
tracking system covers 92% of the solid angle in the cm. 
frame. The track finding efficiency is on average (98±1)% 
for momenta above 0.2 GeV/c and polar angles greater 
than 500 mrad. Photons are detected by a CsI(Tl) elec- 
tromagnetic calorimeter (EMC), which provides excellent 
angular and energy resolution with high efficiency for en- 
ergies above 20 MeV |l^. The energy resolution of the 
EMC is 3% and the angular resolution is 4 mrad for pho- 
tons of energy 1 GeV. 

The asymmetric beam configuration in the laboratory 
frame provides a boost to the T{AS) increasing the mo- 
mentum range of the B meson decay products up to 
4.3 GeV/c. Charged particle identification is provided by 
the average energy loss {dE/dx) in the tracking devices 
and by an internally reflecting ring imaging Cherenkov 
detector (DIRC) covering the central region. A K-n 
separation of better than four standard deviations {a) is 
achieved for momenta below 3 GeV/c, decreasing to 2.5f7 
at the highest momenta in our final states. Electrons are 
identified by the tracking system and the EMC. 

Hadronic events are selected based on track multiplic- 
ity and event topology. We reconstruct B meson can- 
didates from their charged and neutral decay products. 
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including the intermediate states t]' 



{K*^) or i^±^o {Kf^), 



TT 



^0 













The se- 



TT" '^'y^ and 
lection requirements are identical to those used in the 

-T3l. 



branching fraction measurements 12 , 

Candidate charged tracks are required to originate 
from the interaction point, and to have at least 12 DCH 
hits and a minimum transverse momentum of 0.1 GeV/c. 
Looser criteria are applied to tracks forming candi- 
dates to allow for displaced decay vertices. Kaon tracks 
are distinguished from pion and proton tracks via a like- 
lihood ratio that includes dE/dx information from the 
SVT and DCH, and, for momenta above 0.7 GeV/c, the 
Cherenkov angle and number of photons as measured by 
the DIRC. 

We form K^, cj), K *'^, and p*^ candidates from pairs of 
oppositely charged tracks that form a consistent vertex. 
We further combine a pair of charged tracks with a con- 
sistent vertex and a tt" or 77 candidate to select oj or r/'^^^ 
candidates. The Kg candidates are required to satisfy 
|m(7r+7r~) — m/j-o| < 12 MeV/c^ with the cosine of the 
angle between their reconstructed flight and momentum 
directions greater than 0.995 and the measured proper 
decay time greater than three times its uncertainty. 

We reconstruct 7r° (77) mesons as pairs of photons, each 
with a minimum energy deposition of 30 MeV (100 MeV) 
in the EMC. The typical resolution of the reconstructed 
7r° mass is 7 MeV/c^. A ±15 MeV/c^ interval centered 
on the nominal 7r° mass ll^ is applied to select tt" can- 
didates. We combine a p'^ candidate with a photon of 
energy above 200 MeV to obtain an 77^^ candidate. 

We select </), w, 77', and ij candidates with requirements 
on the invariant masses (in MeV/c^) loose enough to 
retain sidebands for later fitting: 990 < m{K^K^) < 



1050, 735 < m(7^^ 



< 830, 930 < m(?77^^ 



990, 900 < to(p7) < 1000, and 490 < m{-fj) < 600. 
The experimental resolutions in the K* and p invariant 
masses are negligible with respect to their natural widths. 
The Kit invariant mass interval is ±150 MeV/c^ for the 
charged and ±100 MeV/c^ for the neutral K* candidates. 
We require the invariant mass of p candidates to be be- 
tween 500 and 995 MeV/c^. 

The helicity angle 6h of a </>, K*, or u resonance is 
defined as the angle between the direction of one of two 
daughters, or the normal to the co decay plane, and the 
parent B direction in the resonance rest frame. To sup- 
press combinatorial background, we require the cosine of 
the K*^ K^TT^ helicity angle, defined with respect 
to the kaon, to be greater than —0.5. This effectively 
requires the vr'^ momentum to be above 0.35 GeV/c. 

We identify B meson candidates kinematically us- 
ing two nearly independent variables |14[ , the energy- 
substituted mass mES = + Pi ■ Pb) /Ef — pg]^/"^ 
and AE = {EiEs - Pi ■ Ps - s/2)/y/s, where {E,,pi) 



is the initial state four-momentum, obtained from the 
beam momenta, and {Eb,Pb) is the four-momentum of 
the reconstructed B candidate. A quantity that is al- 
most equivalent to tties can be obtained from a kine- 
matic fit of the measured candidate four-momentum in 
the T{AS) frame with its energy constrained to that of the 
beam Q . For signal events AE peaks at zero and mES 
at the B mass. Our initial selection requires mss > 5.2 
GeV/c2 and |A.B| < 0.2 GeV. 

Charmless hadronic modes suffer from large back- 
ground due to random combinations of tracks produced 
in the quark-antiquark continuum (e+e^ qq). This 
background is distinguished by its jet structure as com- 
pared to the spherical decays of the B mesons produced 
in the T{AS) decays. To reject continuum background we 
make use of the angle dr between the thrust axis of the 
B candidate and that of the rest of the tracks and neu- 
tral clusters in the event, calculated in the cm. frame. 
The distribution of cos 9t is sharply peaked near ±1 for 
combinations drawn from jet-like qq pairs, and nearly 
uniform for the isotropic B meson decays. Thus we re- 
quire I cos^tI < 0.9 (0.8 for (j>K*'^). We also construct a 
Fisher discriminant that combines eleven variables [ pT| : 
the polar angles of the B momentum vector and the B- 
candidate thrust axis with respect to the beam axis in the 
T{AS) frame, and the scalar sum of the cm. momenta 
of charged particles and photon (excluding particles from 
the B candidate) entering nine 10° polar angle intervals 
coaxial around the i?-candidate thrust axis. Monte Carlo 
(MC) simulation ||l8| demonstrates that contamination 
from other B decays is negligible. 

We use an unbinned extended maximum likelihood 
(ML) fit to extract signal yields and charge asymmetries 
simultaneously. The extended likelihood for a sample of 
N events is 



£ = cxp - ^ 



N 



n-ik 



i,k 



i,k 



(3) 



where Vik [xj ; a) is the probability density function 
(PDF) for measured variables Xj of an event j in cat- 
egory i and flavor state fc, and riik are the yields ex- 
tracted from the fit. The fixed parameters a describe 
the expected distributions of measured variables in each 
category and flavor state. The PDFs are non-zero only 
for the correct final state flavor (fc = 1 for i? — > / and 
k = 2 for B ^ f). In the simplest case, there are two 
categories, signal and background (i — 1,2). The de- 
cays with the charged primary daughter B^ X^^h^ 
{h^ — or TT^, and = 77', w, or 0) are fit simul- 
taneously with two signal (« = 1 for B^ — + X^K^ and 
7 = 2 for B^ X^n^) and two corresponding back- 
ground {i = 3,4) categories. 

We rewrite the event yields riik in each category in 
terms of the asymmetry Ai and the total event yield Ui'. 
mi = X (1 + Ai)/2 and ni2 = x (1 — Ai)/2. The 
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event yields rii and asymmetries Ai in each category are 
obtained by maximizing C ||l^ . The dependence of C on 
a fit parameter rii or Ai is obtained with the other fit pa- 
rameters floating. We quote statistical errors correspond- 
ing to unit changes in the quantity = — 21n(>C/£i„ax), 
where iZmax is the maximum value of the likelihood. The 
90% confidence level (C.L.) limits correspond to a change 
in the of 2.69. When more than one channel is mea- 
sured for the same primary B decay, the channels are 
combined by adding their distributions. 

The PDF Vik {xj ; a) for a given event j is the product 
of PDFs in each of the independent fit input variables 
Xj. These are A.E, m-Es, invariant masses of interme- 
diate states (77', Lu, (j), K*, and 77), Fisher discriminant, 
and the (f> and uj helicity angles for pseudoscalar-vector 
decays. For the simultaneous fit to the decays with the 
charged primary daughter {B^ X^K^ and X^tt^) 
we include normalized residuals derived from the differ- 
ence between measured and expected DIRC Cherenkov 
angles for the . Additional separation between the two 
final states is provided by AE. The separation depends 
on the momentum of the charged primary daughter in 
the laboratory and is about 45 MeV on average varying 
from about 30 MeV for the highest momentum to about 
80 MeV for the lowest momentum primary daughters in 
our final states. 

For the parameterization of the PDFs for AE, m-^s, 
and resonance masses we employ Gaussian and Breit- 
Wigner functions to describe the signal distributions. For 
the background we use low-degree polynomials or, in the 
case of Toes , an empirical phase-space function |Q . The 
background parameterizations for resonance masses also 
include a resonant component to account for resonance 
production in the continuum. In the B decays to vector- 
vector states, the helicity angle distribution is the result 
of an a priori unknown superposition of transverse and 
longitudinal polarizations, and thus is not used for back- 
ground suppression in the fit. For pseudoscalar-vector B 
decay modes, angular momentum conservation results in 
a cos^ 9h distribution for signal. The background shape 
is again separated into contributions from combinatoric 
background and from real mesons, both fit by nearly 
constant low-degree polynomials. The Cherenkov angle 
residual PDFs are Gaussian for both the pion and kaon 
distributions. The Fisher discriminant is described by an 
asymmetric Gaussian for both signal and background. 

The fixed parameters a describing the PDFs are ex- 
tracted from signal and background distributions from 
MC simulation, on-resonance AE and toes sidebands, 
and off-resonance data. The MC resolutions in AE and 
Toes are adjusted by comparisons of data and simulation 
in abundant calibration channels with similar kinematics 
and topology, such as B ^ Dtt, Dp with D Kn, Ktttt. 
The resolutions in the invariant masses of intermedi- 
ate states are obtained from inclusive particle samples. 
The simulation reproduces the event-shape variable dis- 



TABLE I: Results of the ML fits, including number of signal 
events (risig), their charge asymmetry (Acp), and asymmetry 
90% C.L. limits. All results include systematic errors, which 
are quoted second after statistical errors for risig and Acp- 



Mode 






Acp 


90% C.L. 


'/ 

7/777771 


49.5 
87.6 


tin ±3.7 


—Oil ± n 1 1 ± n 02 

-0.17 ±0.15 ±0.01 
-0.05 ±0.15 ±0.03 


Ln 28 4-0 071 




27.6 


1 ?;?±l.9 


-0.01 t ± 0.03 


[-0.50,±0.46] 




31.4 


^li±2.3 


-0.05 ±0.20 ±0.03 


[-0.37,±0.28] 




4.4 
7.1 


i 2:0 ±0.4 
Hi ±1.2 


-0.43 1 °o\fo ± 0.06 
-0.55 1 {5^5 ± 0.05 

r, 01 + 0.54 + 0.10 
U.Oi _ 0.43 - 0.06 


[-0.88,±0.18] 




20.8 


ti:?±i.3 


0.00 ±0.27 ±0.03 


[-0.44,±0.44] 



tributions found in data. The Cherenkov angle resid- 
ual parameterizations are determined from samples of 

£) ° K^TT^ originating from D*^ decays. 

The results of our ML fit analyses are summarized in 
Table |. The signal yields along with branching fraction 
results have been reported earlier [|l^, In all cases 
we find signal event yields with significances, including 
systematic uncertainties, of greater than four standard 
deviations, and hence proceed with asymmetry measure- 
ments. The measured likelihood values are well repro- 
duced with generated samples. The dependence of the 
on Ac p for each decay mode is shown in Fig. 1^ and 
asymmetry measurements are summarized in Fig. g. We 
see no significant asymmetries and determine 90% C.L. 
intervals. 

In the charge asymmetry measurements, systematic 
uncertainties relevant to branching fraction measure- 
ments tend to cancel, but some level of bias is inevitable 
as neither the BABAR detector nor PEP-II is perfectly 
charge symmetric. However these effects are mostly very 
small for the final states considered here. Charge bi- 
ases in track reconstruction and particle identification 
efficiency have been studied in a sample of more than a 
billion charged tracks in multi-hadron events. After pro- 
ton and electron rejection we find an asymmetry in track 
reconstruction efficiency consistent with zero with an un- 
certainty of less than 0.01 for a wide range of momenta 
for tracks originating from the interaction point. Taking 
into account particle identification requirements similar 
to the ones applied to the K* daughters, this consistency 
is still better than 0.02. A D*^ control sample of kaon 
and pion tracks is used to estimate systematic uncer- 
tainties in the asymmetries arising from possible charge 
biases in the Cherenkov angle residual, which are found 
to be less than 0.01. 

From these studies we assign a systematic uncertainty 
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FIG. 1: Dependence of = — 21n(£/£max) on Acp for each 
of the B decay modes. 



of 0.01 on Acp for all the modes with a charged primary 



daughter: rj' 



and (j)K^. For the modes 



with a K* we account for the broader momentum spec- 
trum of the charged daughters and particle identification 
applied to the kaon candidates with a 0.02 systematic er- 
ror. All measured background asymmetries in data and 
signal asymmetries in MC are consistent with zero within 
statistical uncertainties. 

A different type of uncertainty originates in the ML 
fit from assumptions about the signal and background 
distributions. In order to derive systematic errors in the 
event yield and its asymmetry, we vary the PDF param- 
eters with their respective uncertainties. The systematic 
errors in the asymmetries are found to be 0.02 for rj'K^ 



0.03 for W7r± and (t)K^, and 0.06 for (PK*^ 



and 

the latter being dominated by the mode with a 7r°. These 
systematic errors are conservatively estimated and can be 
improved with a larger data sample. 

Uncorrelated (due to PDF variations) and correlated 
(due to selection requirements) systematic errors are 
treated separately in the case of multiple decay channels 
and each is convolved with the likelihood distributions to 
account for all systematic effects in the result. The asym- 
metry measurement in the B ^ ~* (f) K *^ decay mode is 
corrected by the inverse dilution factor 1/(1 — 2ii;), where 
w, the fraction of doubly misidcntified Kt: combinations 
originating from K*'^, is ^0.01. The uncertainties in 
the final results presented in Table ^ are dominated by 
statistical errors. 

In summary, we have searched for direct CP violation 
in charmless hadronic B decays observed in the BABAR 
data. The measured charge asymmetries of the B decays 



FIG. 2: Results of Acp measurements for the B decay modes 
presented in this paper (circles) with the BABAR measure- 
ments in nK modes (squares) shown for comparison 11 
The data sample for the n^K^ result is about 50% larger [11 
Hatched regions are excluded at 90% C.L. 



into final states 77'^=^, t^7r±, (t)K*^, and cp K *° 

are summarized in Table | and Fig. |^. These results, 
along with the asymmetry measurements in i? — > nK 
modes [ p^ 0| and in combination with the earlier mea- 
surements P], rule out a significant part of the physi- 
cal Acp region, allowing for constraints on new physics 
models , but are not yet of sufhcient precision to allow 
precise comparison with Standard Model predictions . 
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